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Abstract: Wc analyse a novel squeezing and entangling mechanism 
which is due to correlated Stokes and anti-Stokes photon forward 
scattering in a multi-level atom vapour. Following the proposal we 
present an experimental demonstration of 3.5 dB pulsed frequency 
nondegenerate squeezed (quadrature entangled) state of light using 
room temperature caesium vapour. The source is very robust and 
requires only a few milliwatts of laser power. The squeezed state is 
generated in the same spatial mode as the local oscillator and in a 
single temporal mode. The two entangled modes are separated by 
twice the Zeeman frequency of the vapour which can be widely tuned. 
The narrow-band squeezed light generated near an atomic resonance 
can be directly used for atom-based quantum information protocols. 
Its single temporal mode characteristics make it a promising resource 
for quantum information processing. 
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1. Introduction 

Quadrature entangled states of light occupy an important place in modern quantum 
information processing, with examples inluding quantum cryptography [I] , teleportation 
[2], computation [3] and last but not least, improving Quantum Non-Demolition (QND) 
measurements [4]. Although continues- variable entanglement and squeezing is a fragile 
resource easily washed out by optical losses, proposals for continuous- variable error 
correction [51 [5] , entanglement purification [7] and distillation [5] have been developed. 

A sub-threshold optical parametric amplifier (OPA) first used to generate quadra- 
ture entangled light in [5] remains the mainstream approach [101 111] . This approach 
offers design flexibility but at the same time requires relatively complicated and sensi- 
tive setup consisting of two or more stabilized cavities. The other possibilities include 
the usage of four wave mixing (4WM) |12| . This and similar methods however generate 
multimode temporal and/or spatial entanglement, that is the output light of the ampli- 
fier contains pairs of frequencies which are independently quadrature entangled. This 
may have adverse consequences for some applications. For example, if multiple modes 
are used in quantum cryptography, modes not detected by the legitimate addressee of 
the information can carry the same signal. In principle they can be detected unnoticed, 
compromising the safety of the protocol. Also frequently a multimode character of the 
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Fig. 1. Off-resonant double A interaction in Caesium in the presence of the mag- 
netic field. The interaction is driven by a strong, linearly polarised probe beam 
detuned by A = 0.85GHz from the extreme atomic resonance. The ground state 
levels are split by Q.l = 322kHz. The collective scattering leads to weak atomic 
excitation to the nif = 3 state described by an annihilation operator b and Q.i 
sidebands of the probe light described by annihilation operators a+ and a- . 

squeezed light is a disadvantage for the protocols involving photon counting, such as 
entanglement purification [5] or preparation of squeezed single-photon states [T51IT^[T5] . 

In this paper we report the proposal and experimental observation of a high purity 
entangled state which at the same time is a pure two-mode squeezed state. The source 
is based on an interaction of off-resonant driving light with the room-temperature spin 
polarised caesium vapour placed in dc magnetic field in relaxation-protected environ- 
ment. The squeezed field is naturally compatible with atomic memories based on the 
same alkali atom as the source. The non-classical state of light is generated in a pair 
of frequency sidebands shifted by the Zeeman frequency around the driving field (the 
carrier) with a unique temporal envelope and the same spatial mode as the driving 
light beam. The frequency of the entangled modes can be tuned by magnetic field and 
the bandwidth of squeezing can be adjusted by changing the driving light power or the 
optical depth of the atomic sample. In the present experiment an ultra-narrow band 
squeezing with approximately 1 kHz bandwidth has been generated. These features 
make it an attractive alternative to existing sources for the applications in quantum 
information processing (QIP), especially those vulnerable to the spurious modes. 

Our method requires almost no alignment and three diode lasers, each of a few mW 
power, two for optical pumping and one for off-resonant driving. 

2. Model of the interaction 

Let us consider an ensemble of atoms (Fig. [1]) with the total ground state angular 
momentum F > 1/2 driven by a strong off-resonant linearly-polarized probe light. Prior 
to the interaction we optically pump all atoms into the extreme magnetic sublevel 
irif = F in the electronic ground state. The orientation is thus orthogonal to both the 
propagation and polarization direction of the driving field. The probe light is seen as a 
superposition of (7+ and G polarizations in the quantisation basis. 

Qualitatively generation of squeezed light by such system can be explained via cor- 
relations between forward scattered TT-polarized photons. The driving light leads to the 
anti-Stokes scattering into mp = F — \ followed by Stokes scattering back to niF = F . 
In the alkali atoms, for the polarization of the driving field perpendicular to the ori- 
entation, the second transition is more probable and therefore almost each anti-Stokes 
photon is accompanied by the Stokes twin with anticorrelated phase. This property is 
mediated by the collective atomic excitation as in the quantum repeater utilizing DLCZ 
protocol [16|. The correlation persists over a period of several milliseconds, the typi- 
cal life tinre of the atonric ground state coherences To in our experiment. If the entire 
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process happens on a timescale shorter than T2, then the emitted Hght is in a single 
mode. The more photon-atom pairs are created and the more atomic excitations are 
released into light, the higher is the degree of squeezing which is defined by the number 
of scattered photons. A general theory of off-resonant interaction of this type neglect- 
ing atomic motion has been discussed in |17| and broadband squeezing was predicted 
in [H]. 

The interaction between light and atoms in our system can be described using two 
creation operators — for the scattered photons and o' for the collective atomic 
excitation into the magnetic sublevel nip = F — I. The Hamiltonian of interaction can 
be written as: 

^int = -tiXaO^b'^ + filpo^b + H.c. , (1) 
where Xa and Xp are the coupling constants, describing Stokes and anti-Stokes scatte- 
ring. At the elementary interaction level XaO^b^ + H.c. describes the active part of the 
interaction, that is photon-atom entanglement, while Xp«^^ + H.c. describes the passive 
part of the interaction, that is a beam splitter - like exchange of excitations between 
photons and atoms. 

In case of realistic multilevel atoms the passive and active coupling Xp and Xa ii^ 
Eq. ([T]) can have substantially different magnitudes. The interaction Hamiltonian can 
be written as Z/int = fiXiPaPb + ^^Xah) where X = Xp + Xa. = {Xp - Ia)l{Xp + Xa) 
14fl2/c'i, fli, fl2 arc the vector and tensor parts of the atomic polarizability, and 
Pa = {d — d^)/iV2 and ph = {b — b^)/iV2 being p quadrature operators associated with 
light and atoms (for details of the Hamiltonian derivation see the Appendix) . For short 
interaction times for alkali atoms with F > 1/2 the second term in Hi^t can be ne- 
glected and the Hamiltonian reduces to the quantum nondemolition (QND) Hamil- 
tonian Hqnd = '^TiXpapb extensively used for spin squeezing, quantum memory and 
teleportation protocols [H [12] . For atoms with F = \/2 the passive and active couplings 
in Eq. ^ have equal magnitudes Xa = Xp and the Hamiltonian is always of the QND 
type. 

However for the stronger coupling/longer interaction time with alkali atoms the com- 
plete Hamiltonian Eq. ([T]) has to be applied which leads to new attractive dynamics. 
If the passive part of the interaction prevails Xp > Z«i which is the case analyzed in 
this paper, £, is real and the interaction leads to swapping of the quantum states be- 
tween light and atoms as discussed below. On the contrary the case of Xa > Xp. i-C- of 
imaginary £, leads to entanglement between light and atoms. 

For a rigorous derivation of the atom-photon dynamics we start with the Hamiltonian 
in the presence of the magnetic field. As illustrated in Fig. [1] the ground state sublevels 
experience Zeeman splitting with the Larmor frequency D.^. Therefore the Stokes scat- 
tering produces a photon in an upper sideband of the probe, described by the creation 
operator a^, while the anti-Stokes scattering couples to a lower sideband of the probe, 
described by the creation operator a[_ . Thus the interaction Hamiltonian is reexpressed 
in the following way: 

^int = -^J^ dz (^-Xadlb"^ + Xpdlb) + H.c. (2) 

L is the length of the cell and we have omitted the space and time dependence of 
operators for brevity. The operator b annihilates an atomic excitation in an atomic slice 
around a given z: 

b= \ -^Y.\mF=^)k{mF^3\k (3) 
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where Na is the total number of atoms in the cell, n is the number of atoms in the 
slice and k indexes the atoms. In the limit of almost all atoms residing in the mp = 4 
state £ is a bosonic operator. = — z'). Analogously. a+ and a_ denote an 

annihilation operator for the upper or lower sideband photonic mode around a given z 
at a certain time with [a-\-(z,t),a_j^{z^t')\ = 5(r — ?'). 

We derive the input-output relations from this Hamiltonian under two assumptions. 
We assume that the light passes through the cell on the timescale much shorter than 
the total time of the interaction or the atomic transient time. Next we assume that 
the atoms see an average over the volume polarization of the light since they move fast 
compared to the evolution of their internal state. Therefore the collective atomic mode 
interacting with light is uniform across the cell with the corresponding annihilation 

operator b = Jq dzb{z) / Vl. We then obtain for the light modes: 

a'_{t)^a_{t) + iXabHt), (4) 

where d±{t) = a±(z = 0,t) and a'j^{t) ~ d±{z = L,t) denote the operators at the in- 
put/output plane of the cell, and we assume that the passage is instantaneous. Note, 
that the state of light inside the cell changes linearly with the coordinate z from d±{t) 
to a±(f)- 

To calculate the atomic state evolution we integrate the Heisenberg equation for the 
local atomic polarisation b over the length of the cell to obtain an equation describing 
the evolution of the entire ensemble: 

= -iXpa+ (0 + iXaai (0 - JswHt) , (5) 

where the factor = |XpP/2— |j„p/2 takes into account the linear change of and 
a_ along the cell, as follows from Eq. ^ . We will see that /sw is the rate at which the 
initial state of the atoms decays and is replaced by the state of the incoming light, hence 
we refer to it as a swap rate. This equation can be readily solved and yields the single 
cell input-output relations. However the result is complicated and not useful for our 
current purpose, because we have disarrayed the simple structure of the Hamiltonian 
by applying the magnetic field. Fortunately, similarly to the QND case this can be 
rectified by letting the light interact with two atomic ensembles in series [501 HI] ■ The 
ensembles are placed in equal magnetic fields, with the atoms pumped into opposite 
extreme magnetic states. Then the upper and lower sidebands interchange their roles. 

To obtain the combined two-cell input-output relations we write the Eqs. ([4]) and ^ 
for both cells, using the output light of the first cell as the input for the second cell. 
Instead of sideband operators fl+ and fl_ we shall use the sine and cosine combinations 
dc = {d+ + d^)/V2 and ds = i{^d++d-) / \/2. For the atoms we use collective annihilation 
operators defined as be = {b\ +b2)/^ and bs = i{—bi +1)2)1^/1 where b\ and ^2 are 
annihilation operators for cell 1 and 2 respectively. We arrive at: 

a''(0 = ac{t) - a/A(f ) + iX.h\{t). 

= - iXpac it ) + iXaa\ (t ) - 7sw^c {t), (6) 

where " denotes the output face of the second cell. We also get an identical pair of 
equations for d^ and bs operators. Below we shall focus on the equations for c operators 
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Fig. 2. Pictorial representation of the interaction. The light in the exponentially 
rising mode at the input of the cell is mapped onto the atoms, and it assumes a 
form of an exponentially falling mode onto which the atoms have contributed. 



implicitly using the fact, that the other solutions will be identical. Integrating the above 
equations gives: 

a^it) =fl,(f) - 2y,^j'^ dt'e-r^-^'-''^a,it') + ie'r^-' (-ZA(O) + zA'(O)) , (7a) 

kit) =ij'^dt'e-y^-^'-''^ {-Zpac{t')+Zaal{t')) +e-1'-'£,(0). (7b) 

Those equations can be cast into a simple form once suitable temporal modes for the 
light are introduced. For probe pulses of duration T we choose: 

fl,(f ) = ^L^Mn.^="' + . . . , m = ^I^A^,„t.-r.w^ + . _ (8) 

where Xl and Pl are quadrature operators for the new modes of light, Mn and A^out 
are normalization factors while the dots stand for contribution from modes orthogonal 
to Nine^''^' and Noute^^"'"^ respectively. We also find it convenient to introduce the 
quadrature operators Xa and Pa for the atomic state: 

With these notations the input-output relations reduce to: 



X'a =^1-^^k^Xa + kPl, Xi ^^T^^Xl + kPa, 



P'a=^/Y~^^Pa-^"kXl, Pi=Vl-^^K^PL-^^KXA. (10) 



where K = 1 ^1 — exp(— 27swr) is the coupling constant. These input-output relations 
have a simple and interesting form. They describe a beam splitter transformation with 
transmission 1 — ^^fc^ and squeezing of the inputs by a factor of 1/^. This is schemat- 
ically depicted in Fig. [21 In particular after a long enough interaction time T S> Ysw 
the initial state of the atoms is squeezed and mapped onto the light, and vice versa. If 
the atoms are initially in the coherent spin state (CSS), the output light will be in a 
squeezed state with a squeezing factor equal 1/^, i.e. Ya,r{P[) =t,^. Note that the above 
equations describe either the cosine combination of the light sidebands and the sym- 
metric excitation of the two ensembles or the sine combination and the antisymmetric 
excitation. Therefore both the P^^. and P[ ^ will be squeezed after the interaction. This 
is equivalent to generating quadrature entanglement between upper and lower sideband 
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Fig. 3. Experimental setup. Cells 1 and 2 are placed in magnetic shields in equal 
magnetic fields B oriented along x. The atoms inside are pumped with circularly 
polarized pump beams of opposite helicity. The probe beam propagates perpendic- 
ular to the B field, interacts with the atoms and finally hits a balanced homodyne 
detector. 



of the probe pulse since 2Var(P£J +2Var(P^ J = Var(P+ + Var(X+ < 2 where 

+ and - denote upper and lower sidebands respectively and the last inequality is exactly 
the quadrature entanglement criterion |24| . 

When tensor polarizability effects can be neglected, Xp — Xa a-nd ^ = 0, the input- 
output relations reduce to the QND case |2l] and entanglement driven by the swap 
interaction disappears. 

The duration of the light mode which contains squeezing can be manipulated by 
changing the driving power and/or the optical depth of our sample. The output mode 
can be also shaped by varying the driving field intensity during the interaction, so 
that 7sw becomes time dependent, which is of particular importance for applications 
for atomic memories. A closer examination of Eq. ((7a|) reveals that the output mode in 
this case has a mode function u{t) proportional to u{t) °^ 7sw(06xp(— /oc/f'7sw(f'))- 
particular, one can shape the driving field pulse such that squeezing is produced in a 
flat top temporal mode. 

3. Experiment 

The experimental setup is depicted schematically in Fig. [31 We use two 22 mm long 
cubic paraffin-coated cells containing about 3.6 x 10^' caesium atoms each. Both cells 
are placed inside magnetic shields in uniform magnetic fields oriented along x axis. 
Each measurement cycle begins with optical pumping of the atoms in the two cells into 
oppositely oriented CSS, as described in [521 [131 [TH] . Next a 15 ms long driving pulse 
blue detuned by 855 MHz from the 6S1/2, F=4^6Pi/2, F'=5 transition is turned on. 
Prior to the interaction it is spatially shaped into a circular flat top beam 20 mm in 
diameter using a telescope beam shaper in order to make the coupling strength uniform 
across the beam. The driving field polarized along the y axis passes through both cells 
along the z-direction. Finally the beam goes through half- or quarter wave plates, a 
polarizing beam splitter (PBS) and onto a pair of balanced detectors. Depending on 
the settings of the wave plates we can either measure the Xi or the quadrature of 
the light. The signals from the detectors are subtracted, sent to a lock-in amplifier and 
digitized with an integrating A/D converter at a 12.5 kHz rate. The number of pumped 
atoms in the cells is monitored by measuring the Faraday rotation of very weak auxiliary 
probe beams propagating through the cells along the magnetic field. 

As a first test, we check the coupling of the mean value of the atomic spin operators 
to the light. The initial spin polarized state (CSS) of atoms corresponds to the vacuum 
state in (Xi, Pa) phase state. A 150 jls long RF-pulse applied orthogonally to the dc 
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Fig. 4. The decay of the mean values {xi{t)) (solid) and {piit)) (dash-dot line) 
measured by the homodyne detector. Prior to this measurement the atoms are 
prepared in a CSS state and subsequently displaced in both and by an RF 
pulse. 

magnetic field in one of the cells displaces the atomic spins equally in Xa and Pa creating 
several tens coherent atomic excitations. After that we detect first Xi and then Pl in 
two series of 200 experimental cycles. From Eq. ()7a|) we can find the expected mean 
values of the output light operators assuming they have zero mean values at the input: 



(4(,)) = v^,-w^P^(, = 0)), 

= -^2^ie-y^-'{XA{t^Q))- (11) 

where the time-dependent light quadrature operators x^(f)-f//7^(f) = ^/2a'l.(t) are written 
with small letters to distinguish them from operators associated with exponential modes. 
The experimental data showing the exponential decay of both mean values as predicted 
by Eq. ([TT|) is presented in Fig.[31 The ratio {xl{Q)) /{pl{0)) found from the figure yields 
(^~^ = 6.3. It agrees very well with the theoretical value at our detuning, which itself is 
not very sharp due to the Doppler broadening. Notice that for a pure QND interaction, 
{piit)) would be zero independently of either measurement time or input mean values 
of the atomic operators. 

Next we proceed to the demonstration of the squeezing of the output light. We prepare 
the atoms in the CSS, send a 15 ms probe pulse and measure the noise in pL,c{t) 
and pL.s{t)i averaging over typically 10000 cycles. Squeezing can be seen in the power 
spectrum of the signal (i.e. {p\ ^{Q.) + p\ ^{£1))) from the lock-in amplifier as shown in 
Fig. [5l When the magnetic field is shifted such that the atomic contribution to the noise 
lies outside the detection bandwidth, we measure the shot noise level, with the spectral 
shape corresponding to the gain/sensitivity function of our detection system. When the 
magnetic fields in both cells are adjusted so that the atoms preccss at exactly 322 kHz 
(the center frequency of the detection range) , we see an apparent dip in the noise power 
at that frequency and around it in the bandwidth of a few hundred hertz. 
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Fig. 5. Average power spectrum {pj^_^{Q.) + pi_,,{Q.)) = Var(/5_|_n + /3_n) + Var(i+£2 — 
of the homodyne signal PL(t). Frequencies around fl^ = 322kHz are shown. 
The atoms are initialised to the CSS state and then a 15 ms long probe pulse is 
shined through them. Solid curve was taken with atoms in both cells tuned to the 
322 kHz Larmor frequency, while dash-dot curve is the shot noise level reference, 
taken with atoms detuned far away from the detection bandwidth. The Gaussian 
shape of the reference spectrum is due to the detection bandwidth, mainly limited 
by the lock-in amplifier used to demodulate the homodyne signal, while the dip 
in the middle is the fingerprint of the ultranarrowband squeezing produced in the 
interaction with the atoms. 

4. Analysis of the temporal modes of the squeezing 

The probe light emerging from the cells has a strong vertically polarised component 
accompanied by a squeezed state in a horizontally polarised component. The latter is 
typically contained in a temporal exponentially falling mode. Since all the squeezing is 
contained in a narrow frequency bandwidth readily detectable by our homodyne setup, 
it is possible to find the excited modes directly from the experimental data. In each 
experimental cycle we sample the output of the homodyne detector at a rate much 
higher than necessary to capture the bandwidth in which the squeezing occurs. This 
way we obtain values of p quadratures of light as a function of time, PL.c{f) Stiid pi.sit)- 
We focus on the cosine quadrature with the results for the sine quadrature being very 
similar. The two-time correlation function {pL,c{t)pL,c{t')) yields the amount of noise in 
any temporal mode or a correlation between pairs of modes that are contained within the 
detection bandwidth. Therefore it is natural to ask whether any mode basis is favoured 
in such case. The answer is provided by the Karhunen-Locve theorem as detailed in the 
Appendix. For each correlation matrix one can find a unique basis of modes M„(f) which 
have no cross-correlations. Each of them can be in principle measured separately and 
it will exhibit a variance of t,^ ~ Var(P^"j) that is also found from the {pL.c{t)pL,c{t')) ■ 
In Fig. [6] we plot the sum of sine and cosine quadrature variances 2Var(P*"j) + 

2Var(P^"j) = Var(P_["^-l-p!"V2 + Var(x|"^-xl"V2 obtained from the measured cor- 
relation matrices. It is expressed as a noise reduction below shot noise. We can see 
that there is one dominant squeezed mode, in which the noise is reduced by 3.5(l)dB. 
Superscript («) denotes n-th characteristic mode function «„(/). Several other modes 
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Fig. 6. Noise reduction lOlogio ^2((P^",')^ + (^^"j)")^ in the characteristic output 
modes. 




Fig. 7. Characteristic output mode functions for the most squeezed mode (blue 
sohd line) , the next mode (green dashed line) and an exponential fit for the most 
squeezed mode (red dash-dot line). 



are squeezed by about 1 dB. This is due to atomic decoherence taking place during 
the interaction. The equations ([7]) describe a perfectly coherent evolution and predict 
pure single mode squeezing. However, in reality the transverse decoherence time T2 is 
comparable to the duration of the fundamental mode. Thus, figuratively speaking, while 
the squeezed light leaves the cells, the state of atoms in the cells is driven back to the 
initial CSS state. This starts the process of squeezing of the light anew in an incoherent 
manner. This is confirmed by the shape of the characteristic mode functions obtained 
from the experiment shown in Fig. [71 The most squeezed mode has an exponentially 
decaying shape, which agrees with the theoretical model. The decay rate is found to 
be (5.5ms)^\ which is virtually equal to the decay rate of the mean values measured 
to be 7sw = (5.7ms)^'. The next mode is rising, supporting the explanation that as the 
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atomic state is brought back to the CSS, an independent squeezing process starts. 
5. Conclusions 

We have developed a model of an ofl-resonant interaction between light and spin po- 
larized multi-level atomic vapour. We find, that for long interaction times and/or small 
detunings, the interaction leads to the swap of states between the light and atoms ac- 
companied by two-mode squeezing (entanglement) transformation of the two sidebands 
of light and the two cells. In particular if the two atomic ensembles arc in a CSS prior 
to the interaction, this state will be mapped with entanglement onto the state of the 
Q.L sidebands of the outgoing light. This state of light is emitted in an exponentially 
decaying temporal mode. At the same time a portion of incoming light, which comes in 
an exponentially rising temporal mode is mapped onto the atoms. 

We have confirmed experimentally that a room temperature Cs vapour generates a 
single temporal light mode in the quadrature entangled state with 3.5(l)dB of entan- 
glement. As predicted by theory the temporal mode in which we find entanglement is 
decaying exponentially. Losses of light and decoherencc of the atomic state populate 
a few other uncorrelated modes weakly squeezed by about 1 dB. We find experimen- 
tally that the rate of the leading dccoherence process scales with optical power and the 
density of atoms approximately in the same way as the rate of the coherent interac- 
tion 7sw. We attribute this decay which currently limits the degree of entanglement to 
light-induced collisions. Applying antireflection coating on the cell windows to reduce 
the losses and reducing atomic dccoherence due to collisions and magnetic dephasing 
should allow to generate even more pure entangled state with a higher degree of squeez- 
ing. We expect that the single mode squeezing and two mode quadrature entanglement 
can be beneficial in QIP protocols where security is a concern. We also expect that it 
can be very useful in protocols where discrimination of a single mode at the detection 
stage is not possible. Note, that the squeezing produced is readily compatible with the 
atomic memories |22] and the mode in which it is produced can be shaped by changing 
the intensity of the local oscillator during the interaction. 

This research was funded by EU grants COMPAS, QAP, and HIDEAS. CM. ac- 
knowledges support from the Elite Network of Bavaria, QCCC. Paraffin coating of the 
cells was skillfully performed by M. Balabas. 

A. Details of the derivation 

It is possible to derive the Hamiltonian from Eq. ^ by evaluating the Clcbsh-Gordon 
coefficients and adiabatically eliminating excited levels. Without further approximations 
this leads to a Hamiltonian describing the scalar, vectorial and tensorial polarizability 
of the atoms, with coupling constants of gq, ai and 02 respectively |19) : 



Above z is the direction of the propagation of the light beam, A is the cross section of 
the beam, 7 is the natural linewidth while A is the detuning. The Hamiltonian p2p can 
be simplified if the atoms are in a state very close to the CSS oriented along x axis. 
The result is different depending on the F number. For concretness we will assume here 
F ~ 4 with almost all atoms in = 4 state with respect to the x-axis, only some in 
mp = 3 and none in the other states. In this situation we can approximate 7^, 7^ and 




(12) 
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7^ by c- numbers and components of j. This way we arrive at: 



Hi, 



'8A 2nAL Jo 



dzl aiS-j, - Ua2{Syjy + 2Sxjx) + ^ 



QQ — a2 \ 16 




The Hamiltonian above may appear comphcated, but in fact only two first terms under 
the integral a — lAajSyjy are nontrivial. The other terms generate classical rotations, 
in particular: the Sxjx term in our settings, when both Sx and jx are macroscopic merely 
causes a rotation in the x-p plane for both light and atoms. The total angle of this 
rotation is typically of the order of a few miliradians and we neglect it. is a constant 
of motion in our system, thus the term jx only shifts the Larmor frequency — it 
represents the Stark shift. The last term, 5^, adds to a rotation of in the x^-pc plane 
for light but is still negligible. 

Finally we can rewrite nontrivial terms from the above Hamiltonian a\Svj~— \Aa2Syjy 
using bosonic operators, a for light and b for atoms. This is accomplished using the 
relations: 



Sy=^{a + a^) 
'4 + b') 



Jz 




(14) 



where 't' is the photon flux per unit time and Na is the total number of atoms and 
the sign in S, is a consequence of negative sign of Sx in our settings. This way we can 
approximate the interaction Hamiltonian in the form: 



I b-b'' 



iV2 



.b + b' d + a' 



(15) 



where 



7sw = 14flifl2 



A2 UA27r 



1402 



(16) 



The expressions for the vector a\ and 02 tensor polarizabilities for Cs D2 line can be 
found, for example in |21j . It can be directly verified, that the above Hamiltonian is 
identical with Eq. © , with Xa = V7sw/2( 1 - ^ ) and Xp = VlW2 ( 1 + ^ ) • 

Let us note that the sign of can be flipped, that is the Xa and %p can be inter- 
changed. This is accomplished by rotating the polarisation of the driving field by 90°, 
Sx —Sx or by switching from blue to red detuning A —A. The input-output relations 
given in Eq. pU]) remain valid for both signs of However for a real ^ they describe 
the entanglement between two sidebands of the light field, whereas for imaginary ^ they 
entail the entanglement between the light and atoms. 



B. Eigenmode decomposition 

In the experiment we measure directly the covariance matrix of the homodyne signal 
Pcit): 

C{t,t')^{p,{t)pc{t'))~{p,{t)){p,{t')) (17) 
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where /?c(f) is a properly scaled signal directly form the lock-in amplifier. After the 
measurement we can calculate the amount of noise in any temporal mode characterised 
by the modefunction M„(f), where n indexes a set of modes we are interested in. A P- 
quadraturc operator for the nth mode is simply Qn = J u„{t)pc{t). We can calculate both 
the variance of any Q„ and correlation between any two of them: 

{QnQ,n) - {Qn){Qm) = j j dtdt' U„{t)C{t ,t')Un,{t') (18) 

According to the Karhunen-Loeve theorem one can find a set of mutually uncorrelated 
modes u„{t) by simply performing a spectral decomposition of measured C{t,t'). In this 
way we obtain eigenvalues ^„ and eigenfunctions u„(t): 

C{t,t')^Y.^nUn(t)Un{t') (19) 
n 

The eigenvalues B,n are equal to the quadrature variances, t,„ = {Q^) — {Qn)^, while the 
eigenfunctions give the quadrature mode functions. The same procedure can be repeated 
for the sine component of the homodyne signal, Ps{t) and it yields results identical to 
within the experimental uncertainties. 
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